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Abstract 
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guarantee the stability of a dark matter candidate and to forbid the flavour changing neutral current at 
the tree level, respectively. The dominant form factor Fz of the H^W^Z vertex can be enhanced by 
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I. INTRODUCTION 


Although the discovery of the Standard Model (SM) like Higgs boson at the Large Hadron 


Collider (LHC) [l 


^ suggests that there is an isospin doublet scalar field in the Higgs sector, the 


possibility of the existence more Higgs doublets still remains open. In fact, a second doublet is 


B- 


often introduced in new physics models such as the Minimal Supersymmetric SM (MSSM) 

In addition, models with a multi-doublet structure have also been discussed based upon various 

n 

physics motivations, e.g., to emlain tiny neutrino masses via radiative generation [g], to provide a 
dark matter (DM) candidate B and to supply extra CP violating phases [s] for the explanation of 
the baryon asymmetry of Universe. Thus, testing the existence of additional doublet fields is quite 
important to probe new physics scenarios beyond the SM. 

One of the most important features of models with multi-Higgs doublets is the appearance of 
physical extra scalar bosons such as charged Higgs bosons In particular, the properties of 

states strongly depend on the structure of the Higgs sector, e.g., the symmetries of the model, 
the actual number of doublets, the mass spectrum, etc. Therefore, through the detection of 
and by measuring those properties, e.g., the mass, couplings, production cross sections and decay 
rates, one can directly probe the existence of additional doublets as well attempt extracting the 
structure of the Higgs sector. 

Among the various observables related to , studying the H^W^Z vertex is quite interesting 
because of the following features. Firstly, it has been known that the H^W^Z vertex does not 


appear at the tree level^ in multi-doublet models B) because of an approximate global SU{2) 
symmetry known as the custodial symmetry^ in the kinetic terms for the doublet fields. Secondly, 
although the H^W^Z vertex is loop induced, its magnitude can be enhanced by non-decoupling 
effects of particles running in the loop, especially for the case where they come from the sector 
which breaks the custodial symmetry. For example, the top and bottom quark loop contributions 
to the H^W^Z vertex give the quadratic dependence upon the top quark mass [l^, which is 
responsible for the violation of the custodial symmetry in the Yukawa sector. In Refs. 0,0], 
the impact of extra Higgs boson loop contributions on the H^W^Z vertex has been evaluated in 


the 2-Higgs Doublet Model (2HDM) B- It been shown that a large mass splitting between 
the CP-odd Higgs boson and the charged one gives a sizable correction to the H^W^Z vertex. 


^ The H^W^'y vertex does also not appear at the tree level in any models with the 17(l)em symmetry. 

^ In fact, the custodial symmetry is broken by the U{1)y coupling in the kinetic sector which generates the mass 
difference between the W and Z bosons. 
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From the above reasons, it is clear that the strength of the Z vertex measures the effects of 

the violation of the custodial symmetry in the model embedding it. Therefore, by measuring this 
vertex, we can indirectly observe such a new physics effect. 

Feasibility studies to measure the H^W^Z vertex have been performed in Ref. [l^ for the LHC 


and in Ref. 


14i | for future linear colliders. 


In this paper, we calculate the magnitude of the H^W^V {V = Z,^) vertices at the one- 
loop level in the 3-Higgs Doublet Model (3HDM), in which the Higgs sector is composed of two 
active (with a non-zero Vacuum Expectation Value (VEV)) and one inert (without a non-zero 
VEV) doublet scalar fields. In this model, the scalar bosons from the inert doublet field give an 
additional contribution to the H^W^V vertex with respect to the top/bottom quarks and scalar 
bosons from the active doublet loop contributions. As a phenomenological application, we also 
discuss how such new contributions change the decay branching fractions of the —>• W^Z 
and 7 modes and, consequently, the production cross sections involving these decay 

processes at the LHC. 

This paper is organized as follows. In Sec. m we define the Lagrangian of the 3HDM, i.e., the 
scalar potential and the Yukawa interactions. In Sec. imi we introduce the form factors of the 
H^W^V vertices and discuss relationships between these form factors and effective operators. We 
then explain how to calculate these form factors at the one-loop level. In Sec. IIVI we summarise 
various constraints on the parameters of our model. From the theoretical point of view, we consider 
vacuum stability and perturbative unitarity. As experimental constraints, we take into account the 
bounds from the Electro-Weak (EW) 5, T and U parameters, the flavour experiments and direct 
searches for states from LEP-H and the LHC Run-I. In Sec. El we show numerical results 
for the form factors of the H^W^V vertices, branching fractions of and their signal cross 
sections at the LHC. Our conclusion is given in Sec. IVIl In Appendix, we present the full analytic 
expressions for the form factors of the H^W^V vertices. 

II. THE MODEL 


We give a brief review of the 3HDM^ of which the Higgs sector is composed of two active and 
one inert isospin doublet scalar fields [l5, 16|. We represent the active doublets as <I'i and $2 
whereas the inert doublet as 7 . Such an inert nature can be realised by assuming an unbroken Z 2 


The model with two inert plus one active doublets have been discussed in Refs. |l7l. Il8||. 
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{Z2, Z2) charge 

Mixing factor 


4 >i 4)2 ry Ql Ll ur (Ir cr 


Type-I 

(+)+) (+)“) (“1+) (+1+) (+1+) (+i“) (+i“) (+i“) 

cot (3 cot P cot P 

Type-II 

(+)+) (+)“) (“:+) (+1+) (+1+) (+:“) (+1+) (+1+) 

cot P — tan P — tan P 

Type-X 

(+)+) (+)“) (“1+) (+1+) (+1+) (+i“) (+i“) (+1+) 

cot P cot P — tan P 

Type-Y 

(+)+) (+)“) (“1+) (+1+) (+1+) (+i“) (+1+) (+i“) 

cot P — tan P cot P 


TABLE I: Charge assignments of the unbroken Zi symmetry and the softly-broken symmetry. The 
mixing factors in the Yukawa interaction terms in Eq. (HI are also shown. 


symmetry in the scalar potential, in which only ry has an odd parity while all the other fields are 
assigned to be even. One of the important consequences of imposing such a Z 2 symmetry is that 
the lightest neutral scalar component in r] can be a DM candidate, because it cannot decay into 
SM particles. 

In addition to the Z 2 symmetry, we impose another Z 2 symmetry, denoted by Z 2 to distinguish 
it from the above one, which is required to forbid the Flavour Changin g N eutral Current (FCNCs) 
at the tree level. This prescription is the same as that in the 2HDM |19l ]. For the Z 2 symmetry, 
we consider the softly-broken case, since avoidance of tree level FCNCs can already be achieved 
in this case. Under the Z 2 symmetry, four independent types of Yukawa interactions (Type-I, -II, 
-X and -Y) 2rH22l| are allowed depending on the assignment of the Z 2 charge to the SM fermions. 
In Tab. HI we show the charge assignments required by the Z 2 and Z 2 symmetries for the three 
scalar doublets $ 1 , $2 and ry and all the SM fermions, where (e^) is the left (right)-handed 
lepton doublet (singlet) and Ql dj^) is the left (right)-handed quark doublet (up-type and 
down-type quark singlets). 
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A. The scalar potential 


The most general scalar potential under the SU{2)l x U{ 1 )y x Z 2 x Z 2 symmetry is given by 


V{<^i,<l> 2 ,v) = + h.c.) 

+ + iA2(^^$2)^ + A3 ($I^>i)($^^ 2) + A4|^>1$2|^ + ^[A5(^>1^'2)^ + h-C.] 

+ ^Kiv^v? + Pii^i^i)iv^v) + P 2 \^\ri\‘^ + + h-c.] 

+ ai{<^l^ 2 ){p^v) + + h.c.], (1) 

where /X 3 , A 5 , ps and 173 are complex parameters in general. Throughout the paper, we take these 
parameters to be real for simplicity. The scalar fields can be parameterised as 


= 


w- 


+ Uj + izi') 


, (i = l,2), r] = 


. 75(^iT + . 


( 2 ) 


where Vi are the VEVs of <hj with ~ (246 GeV)^. The ratio of the two VEVs is 

parameterized as the usual way by tan (3 = V 2 jvi. 

The mass formulae for the active sector are exactly the same as those in the 2HDM at the tree 
level. The mass eigenstates for the active scalar bosons are given as: 



R{9) 


( cos 9 — sin 9 
sin 9 cos 9 



( 3 ) 


where and are the Nambu-Goldstone (NG) bosons which are absorbed as their longitudinal 
components by the and Z bosons, respectively. We define the h state to be the SM-like Higgs 
boson with a mass of about 125 GeV discovered at the LHG. 

The squared masses of the and A states are then calculated as 

y2 

m^± = - y(A4 + As), = As, (4) 


where 


AR = 


pI 


sin f 3 cos /3 


( 5 ) 
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The squared masses for the CP-even scalar states and the mixing angle a are expressed by 


(6) 

( 7 ) 

( 8 ) 


m\j = cos^(a — + sin^(a — /3)M|2 + sin2(a — j 5 )Ml 2 , 

m\ = sin^(a — / 3 )Mii + cos^(q: — /3)M|2 — sin 2 (Q; — / 3 )M^ 2 i 

/ 2M?2 

tan2(Q: — ( 3 ) = 




where Mf- {i,j = 1 , 2 ) are the mass matrix elements in the basis of {h'^, h'2) defined in Eq. 


= u^(Ai cos'^ /3 + A2 sin'^ / 3 ) + -^(Aa + A4 + A5) sin^ 2 / 3 , 


M22 = sin^ /3 cos^ /3 [Ai + A2 — 2(A3 + A4 -|- A5)], 


..2 „■ 2 


Mf2 = — sin 2 / 3 (—Ai cos^ /3 + A2 sim ( 3 ) + — sin 2/3 cos 2/3(A3 + A4 + A5). 


( 9 ) 


Because of the unbroken Z 2 symmetry, the scalar bosons from rj do not mix with those from d>i 
and <f> 2 . Thus, the mass formulae of the inert scalar bosons are simply given by 

^2 

m^± = y [pi cos^ /3 + ai sin^ / 3 ] , ( 10 ) 

^2 

+ — [(pi + p2 + /03) cos^ /3 + (fJi + 0-2 + 0-3) sin^ ( 3 ] , (11) 

^2 

= p,‘^ + — [(pi + p2 - Pa) cos^ /3 + (fJi + 0-2 - 0-3) sin^ ( 3 ] . (12) 


B. The Yukawa Lagrangian 

The most general form under the Z 2 symmetry is given by 

—Cy =YuQ^ dQL^ddji + YeLl^e^R + h.C., (13) 

where ^u,d,e are $1 or $ 2 - The interaction terms are expressed in terms of mass eigenstates of the 
Higgs bosons as 

-4“= E + 

f=u,d,e 

\Pi 

+ — [VudU {mdU Pr - muiuPi) dH^ + me^evPReH^ + h.c.] , (14) 

where 4/ is the third component of the isospin for a fermion /. In Eq. (1141) . and are defined 
by 

= sin(/3 - a) + C/ cos(/3 - a), 

= cos(/3 — a) — sin(/3 — a), 
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( 15 ) 

( 16 ) 



and in each type of Yukawa interactions are listed in Tab. HI 

It is important to mention here that there is the so-called SM-like limit or alignment limit 
defined by sin(/3 — a) 


231, 


24( 1 ■ In this limit, all the h coupling constants to the SM particles 


become the same values as those of the SM values. In fact, the ratios of hff and hVV couplings 
in our model to those in the SM are respectively given as given in Eq. m and sin(/? — a). 

III. THE H^W^V VERTEX 


The amplitude of —)• W^V {V = Z, 7 ) is expressed as 

iM{H^ ^W^V) =igmwV^''ewij,{pw)^Vu{pv), for Y = Z, 7 , 
where Vy'^ is written in terms of the following three dimensionless form factors; 

= g^^’^Fv + Hv, 


mi 


o n nrj^ TVio mm o' t'oOT rb o ■Povt* T/T/'^ 


(17) 


(18) 


with and py being the incoming momenta for and V, respectively. For the case of Y = 7 , 
the Ward identity guarantees the following relation; 


= 0 . 


From this relation, the form factor is written as 


F^ = 


_ G-y 


1 - 




m 


w 


(19) 


( 20 ) 


where we use and (pyr + Pj)^ = m^±. 

In our model, the H^W^V vertices do not appear at the tree level, just like in the 2HDM. This 
is clearly seen by introducing the so-called Higgs basis of the active scalar doublets defined as 


$1 


$2 


= R{P) 




( 21 ) 


where 


d> = 


G+ 

Mh'^+v + iG^) 


L v /2 


T = 


H+ 


-^{h '2 + *^) 


( 22 ) 


with/i'^ = 77 cos(/3 —a)-|-h sin(/3 —a) and /12 = —77 sin(/3 —a)-|-/i cos(/3 —a). The kinetic Lagrangian 
for <I>i and ^2 is then rewritten as 


Gkin = = 177^^1" + 


(23) 
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where is the covariant derivative. Since the gauge-gauge-scalar type vertex is proportional 
to the Higgs VEV v, these vertices come from the term as only $ has a non-zero VEV. 

However, the physical charged Higgs bosons are contained in the term. Therefore, the 

H^W^Z vertex is absent at the tree level^. The above statement can be generalised to a model 
with N active doublet scalar fields. In that case, we can also define a base transformation similar 
to the one of Eq. (EH). Regarding the H^W^'y vertex, it does not appear at tree level in any 
models based on the SU{2)l x U{1)y —>• C/(l)em gauge theory, because of the [/(l)em invariance 
and the consequent Ward identity. 

The form factors dehned in Eq. ()18p are introduced from the following effective Lagrangian [l^ . 
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Aff = fzH+W-Z>^ + gyH+F^Fv^,u + ihye^,p„H+Fi;^ F^^ + h.c., (24) 

where F^ and Fy^ are the field strength tensors for and V, respectively. It can be seen that 
the coefficient fz has mass dimension one whereas Qu^wv and hz have mass dimension minus 
one. Hence, the coefficient fz can be proportional to a squared mass (M?) of a particle running 
in the loop according to a dimensional analysis: 

/\ y'2 

fz--ggz^HMf), (25) 

V 

where is a dimensionless function. Typically, it is expressed by the logarithmic function of Mf. 
On the other hand, gz and hz can be expressed as 

gz, hz ^ (26) 

V 

where Q is another dimensionless function of Mf. Therefore, only the coefficient fu^wz can be 
enhanced significantly due to the Mf dependence, so that the form factor Fz gives the domi¬ 
nant contribution to the H^W^Z vertex. In fact, it has been pointed out in Ref. [3 that the 
top/bottom loop contribution to the form factor Fz is proportional to mf only, as mt 3> m^. The 
origin of the quadratic dependence can be understood in terms of the Yukawa coupling H~^tb, which 
is proportional to mt/v as in Eq. ()14|1 . and of another mt coming from the chirality flipped effect. 
Similarly, the quadratic mass dependence appears in the extra Higgs boson loop contribution as 


* If we consider models which contain scalar fields with isospin larger than 1/2 such as tripletSjthe Z vertex 

can appear at tree level. The expression for the H^W^Z vertex can be found in Refs. [3, Q in the general 
extended Higgs sector which contains Higgs multiplets with the isospin T and the hypercharge Y. In addition, 
it has been known that in models with an extension of the gauge sector such as SU(2) x SU(2) x 1/(1) [2^, the 
H^W^Z vertex also appears at the tree level. 





discussed in Ref. [ll|. This too can be understood, as the trilinear H^SS' {S and S' being extra 
scalar bosons) couplings can be rewritten by squared masses of extra scalar bosons. 

Another important reason for the appearance of a Mf dependence in Fz is in relation to a 
violation of the custodial SU{2)v symmetry. As it has been discussed in Ref. [ll|, the dimension 
three term in Eq. (1241) comes from the following operator® 




(27) 


where $ = (<1>'^, <h) and T' = T) with = ia 2 ^* and = ia 2 ^* are the 2x2 representation 
form of the Higgs doublets. They are translated under the SU{2)i x SU{2)ii symmetry by $ —)• 
and ^ where Ul and Ur are respectively the SU{2)l and SU{2)r unitary 

transformation matrices. We can see that the operator given in Eq. (1271) is not invariant under the 
SU{2)r transformation, so that this operator breaks the SU{2)r invariance. Since the custodial 
SU{2)v symmetry corresponds to the remaining symmetry after the EW symmetry breaking, i.e., 
SU{2)l X SU{2)r — )• SU{2)v and a violation of the SU{2)r symmetry means a violation of the 
SU{2)v symmetry. 

Therefore, the quadratic mass dependence in Fz can be understood as a result of the custodial 
symmetry breaking. In fact, it has been known that the mass difference between the top and 
bottom quarks gives the violation of the custodial symmetry in the Yukawa sector. In addition, that 
between A and also gives the violation of the custodial symmetry in the Higgs potential 26]. 


Since the top quark mass is already known by experiments, the top quark loop contribution to 
the H^W^Z vertex is determined by its mass®. In contrast, parameters in the scalar sector have 
not yet determined by experiments except for the Higgs boson mass of about 125 GeV, so that we 
can expect a sizable enhancement of the H^W^Z vertex from scalar boson loop effects in suitable 
regions of the 3HDM parameter space. 

In the following, we discuss how we calculate the form factors of the H^W^V vertices. We 
can separately consider the one-loop contributions to the vertices from the IPI diagrams and the 
counter terms as 


(Fy, Gy, Hv) = + 5Fv, GiF^ + 5Gy, + 5Hv), (28) 


® The operator also gives the Z term in the effective Lagrangian which is proportional to 

sin^ Ow- However, such an effect is cancelled by the counter term of the H^WZ vertex. 

® In our model, the top quark loop contribution also depends on tan /3, and in all the four types of Yukawa interactions, 
its dependence is given by cot /3. 
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where and SXy are respectively the IPI and the counter term contributions to the form factor 
Xv = F, G and H). Their analytic expressions are given in App. A. 

The counter term contributions are obtained as follows. First, we dehne the renormalized two 
point function for the W^-H^ mixing as 

^wniP^) = i-ipn'^WHip^), (29) 

where is the incoming four momentum of . The renormalised form factor T wh is given by 

^wh{p‘^) = imw^GH + ^wnip'^)^ (30) 

where 6gh is the counter term for the G^-H^ mixing, and T^]^ is the IPI diagram contribution 
to the W^-H^ mixing. The analytic expression of is given in App. A. The counter term is 
obtained by the shift of the charged NG boson held G^\ 


G 


'± 


(1 + 5Zg/2)G^ + 6ghH^ 


(31) 


By imposing the on-shell renormalisation condition [27], [28 


^wh{p^ = mjj±) = 0 , 


(32) 


we can determine the counter term 


6gh = i- 


mw 


We then obtain the counter term contribution to the H^W^V vertex as 


(33) 


P-GWV —- 'iGi\ys‘^G~^W^ em\YG^W h.c. 




- mw s^Sgw n + emw^GwH'^W^ + • • •, 


-w 


(34) 


where = sin^w and = cos^w with 9w being the weak mixing angle. From Eqs. (1331) and 
(IHl), 6Fv is given by 




mw 


mw 


We then obtain the hnite results for the form factors of the H^W^Z and Fl^W^'j vertices. In 
the case of sin(/3 — a) = 1, mj^ = mA S> mf{± and S> m^±, we obtain 


Fz 


cot /3 


IGvr^u^c 


w L 


Ncmj + 


m2 - 


m, 


(tan l3-l) + [m^^- m^± - —p 2 


(36) 
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where the hrst, second and third terms correspond to the contributions from t-b, active and inert 
scalar boson loops, respectively. From the above expression, we can clearly see the quadratic mass 
dependences mf, and However, as it will be discussed in the next section, the case 

considered in the above, i.e., rrijj = ruA mfj± and S> m^± also gives the similar 

quadratic dependence in the EW T parameter. Therefore, too large mass difference between 
and A (with and that between and (with are not allowed. Instead 

of taking the above case, we can consider the case with sin(/3 — a) = 1, tua S> mfj±{= m^) and 
^ ^??±(= where the contribution to the T parameter from extra scalar boson loops is 

cancelled. We then obtain 


Fz ^ 


cot /3 


Ibvr^u^c- 


w 


Ncmj + (M^ — m^±)(tan^ [3 — 1)F 


(^) 

\m\ J 


- -i^{p2 + P^)F 




m. 


VA 


where W = 3 is the color factor, and the function F is given by 

F(r) = -^- — [3 — 4r + + 2(2 — r)r Inrl-Inr. 

^ ^ 4(1 -r)2 L ^ ^ J 2 

This function has the following asymptotic behavior: 

3 1 1 1 — r 

F(r) ~-Inr for r ^ 1, F(r) ~ — for r S> 1, F(r) ~- for r ~ 1. 

42 ^ 4 ^ ' 2 


(37) 


(38) 


(39) 


In this case, although the quadratic dependence m\ and disappears, there still remains their 
logarithmic dependence. 


IV. CONSTRAINTS 


A. Vacuum stability 


The stability condition for the Higgs potential is given by requiring that the Higgs potential 
is bounded from below in any direction of the scalar boson space. The necessary and sufficient 


161 as 


condition to guarantee such a positivity of the potential has been derived in Ref. 

>0, Ai > 0, A2 > 0, (40) 

\/A1A2 + A > 0, AiyAi + p > 0, Y^AjyA^ + h > 0, (41) 

A^A + \/Aict + X 2 P > 0 or A^A^ + Aih^ + X 2 p‘^ ~ A^AiA2 ~ 2Xpa < 0, (42) 

A = A3 + MIN(0, A4 + A5, A4 — A5), 
p = Pi + MIN(0, p2 + P3, P2 - Pa), 

h = (Ti + MIN(0, (T 2 + (T 3 , 0-2 —(T 3 ). (43) 
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B. Unitarity 


Some combinations of scalar quartic couplings are constrained from perturbative unitarity. In 
the 3HDM, the s wave amplitude matrix for all the 2-to-2 body scalar boson elastic scatterings 
have been calculated in Ref. 29| in the high energy limit. We obtain the following independent 


eigenvalues or sub-matrices for the s wave amplitude matrix as 



^ 3 Ajj 2piP 2 2 ( 7 i -|-( 72 '^ 


^ \ P2 0 - 2 ^ 


P3 0-3^ 

Xi = 

2pi p 2 3 Ai 2A3 -|- A4 

, ^2 = 

P 2 Ai A4 

, ^3 = 

P3 Ai As 


y2(7i -|- (72 2A3 -|- A4 3A2 j 


(^(72 A4 \2j 


\ 0'3 As X 2 J 


2 /f = As-I- 2 A 4 ± 3 A 5 , (45) 

yf = Pi + 2p2 ± 3 / 93 , (46) 

= fji-b 2(72 ± 3(73, (47) 

2 /J = A 3 ±A 5 , (48) 

yf = Pi ±P3, (49) 

//^ = cri±(73, (50) 

Vj = A 4 , (51) 

y^ = pi±P 2 , (52) 

yg=ai±a2. (53) 

We then require the following condition: 

Ixil < Svr, \yf\ <871, (i, j = 1, ...9), (54) 

where Xj are the eigenvalues of Xi, X 2 and X^. 


C. S, T and U parameters 

The EW oblique parameters S, T and U Q can be modihed from the SM prediction by 
the extra scalar boson loop contributions and the modified SM-like Higgs boson couplings. The 
differences in the predictions of the S, T and U parameters in the 3HDM and those in the SM are 
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FIG. 1: Constraint from the S and T parameters on the plane in the case of uriH = ruA, sin(/3—a) = 

1 and m^± = = mAl‘2.. The charged Higgs boson mass is fixed to be 150 GeV (left panel) and 200 

GeV (right panel). The 95% CL excluded regions are indicated in the figure. 


given in the case with sin(/3 — a) = 1, as 


AT- 


247r2 


:(m^± -niAf, 


AU 


Utt 


m?A 2mrj± 

In ^ 5 ^ +-^ - 2 ~ 0, 


m 




ruA 


assuming ~ 


IZTT V 


f In + In , for > m^±, 

V m„j_ 6/ ^ 

r/± 


AS ~-(In 


m, 


127r V 


A5~ 


127r 


In ■ 


m, 


5 

6 

m. 


for 


+ 




m 


H± 


The general expression is given in Ref. 
AT are extracted by fixing AU = 0 as 


m I 

T?± 


— 1) ~ 0, for m^± ~ . 


(55) 

(56) 

(57) 

(58) 

(59) 


29l |. From the global fit of the EW precision data, AS and 


AS = 0.05 ± 0.09, AT = 0.08 ± 0.07, 


(60) 


with the correlation coefficient of +0.91 


3ll 


In Fig. [H we show the constraint from the S and T parameters on the plane. We take 

sin(/3 — a) = 1, and = mAl‘2.i which is also taken in the numerical results 

shown in Sec. lYl In the left and right panel, mH± is fixed to be 150 GeV and 200 GeV, respectively. 
We can see that, for ~ , a magnitude of the mass splitting between A and to be larger 

than about 75 GeV is excluded by the T parameter due to the quadratic dependence of the mass 
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splitting shown in Eq. (I55p . In this case, the contribution to AS is almost zero as it is seen in 
Eq. (ISUI) . Conversely, in the case of S> m^±, the positive logarithmic contribution to AS 
appears as shown in Eq. (1571) and a too large mass splitting between r]^ and r]^ is excluded by AS". 
However, the constraint from AS* is getting milder when there is a positive contribution to AT, 
because of the positive correlation between AS* and AT. Therefore, in order to have a large mass 
splitting between and rj^, which is required to obtain a significant contribution to the H^W^Z 
vertex, we need a mass splitting between A and H^. 


D. Flavour constraints 

We can apply the same constraints from the B physics measurements as those in the 2HDM 
to our 3HDM, because of the same structure of the active sector. From the 6 —)• sy process, 
the mass bound of mi{± > 322 GeV is given at 95% confidence level (CL) in models with the 
Type-II and Type-Y Yukawa interactions with tan (5 > 2 via the next-to-next-to-leading order 
calculation performed in Refs. [^, 33]. This bound is getting stronger when a smaller value of 
tan [3 is considered. In models with Type-I and Type-X Yukawa interactions, the constraint from 
6 —?■ S 7 is only important in the small tan/3 case. For instance, the lower limit on mij± is given to 
be about 100, 200 and 800 GeV at 95% CL in the cases of tan/3 = 2.5, 2 and 1, respectively 

The B^-B^ mixing also gives a bound on m}j±, especially for small tan ,0’s. In the case of 
tan/3 = 1, rnu± < 500 GeV is excluded at 95% CL in models with all the types of Yukawa 




interactions 


34( 1 . which is stronger than the constraint from 6 —)• sy for the Type-II and Type-Y 


cases. This bound becomes rapidly weaker when we consider tan/3 > 1, e.g., for tan/3 = 1.5 (2), 
the limit is < 300 (100) GeV at 95% CL. 

E. Direct search at LEP II 


I experiment, charged Higgs bosons have been searched via the e"’'e —)• Z* 


35( 1 . From the non-observation of a significant excess, the lower mass limit has been 


At the LEP 
H'^H~ process 

taken to be about 80 GeV at 95% GL under the assumption of BR(LI^ —)• t^v) + cs) = 

1. The slightly stronger bound mfj± > 90 GeV can be obtained assuming BR{H^ —)• r^z/) = 1. 
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FIG. 2: (Left panel) The product of branching fractions BR(t —>■ iJ+6)xBR(iL+ —r+z/) as a function of 
tan P in the Type-I (red curves) and Type-X (blue curves) 2HDMs/3HDMs. We take TO_f/± = = mn = M 

and sin(/3 — a) = 1 in this plot. The dashed and solid curves respectively show the cases of TOj/i =100 
GeV and 150 GeV. The horizontal dotted lines show the upper limits (0.23% and 1.3%) from the LHC data. 
(Right) Excluded parameter regions on the iaii P-mH± plane in the Type-I and Type-X 2HDMs/3HDMs. 
Regions inside from each curve are excluded at 95% CL by the measurement of top decay t —>■ H^b —>■ r^bv. 
The solid and dashed curves are the results using the upper limit on BR(f —iL+6)xBR(iJ+ — t^v) to be 
0.23% and 1.3%, respectively. 


F. Direct search at LHC Run-I 


At the LHC, searches have been performed for the two cases: the low mass region m[j± < 
nit+mb and the high mass region mjj± > mt+rrib. For the low mass case, the t —)• H^b decay is used 
as the production mode and the full process pp ^ ti ^ hbH^W^ with the decay 

has thus been analysed. Using the data obtained at ^/s = 8 TeV after 19.5 fb~^ of the integrated 
luminosity, the upper limit on the product of branching ratios BR(t —)• Ff^6)xBR(iL^ —)• t^v) 
has been obtained to be between 0.23% and 1.3% at 95% CL for m}j± in the range of 80 GeV to 


160 GeV 


361] 


In the left panel of Fig. [2l the above product of branching ratios is shown as a function of tan /3 
in the Type-I and Type-X 2HDMs. Because the light scenario, i.e., m[j± < nit, in the Type-II 
and Type-Y 2HDMs has already been excluded by 6 —)• sy data as explained in Sec. IV-D, we 
here only show the Type-I and Type-X cases. In the Type-X 2HDM, the product of the branching 
fractions is slightly larger than that in the Type-I 2HDM. This can be understood in such a way 
that in the Type-X 2HDM the branching fraction of —)• is enhanced as tan/3 is increased, 

while it does not depend on tan/3 in the Type-I 2HDM. For example, BR(FI“’" —)• t~^v) can be 
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almost 100% when tan /3 > 3 in the Type-I 2HDM, but it is about 40% in the Type-I 2HDM. In 
contrast, the branching ratio of t —?■ H^b is given by the same value in both Type-I and Type-X 
2HDMs. Therefore, a bit stronger bound on tan/3 for a fixed value of mfj± is obtained in the 
Type-X 2HDM. For example, if we use the stronger bound for BR(t —)• —)• t^u), 

i.e., 0.23%, tan/3 < 6 (4) and 15 (10) are excluded for mfj± = 100 and 150 GeV in the Type-X 
(Type-I) 2HDM. 

For the high mass region, i.e., mfj± > rrit, the production process gb —)• tH^ (i-e., H^- 
strahlung) can be used instead of the top quark decay^. The 95% CL upper limit on the cross 
section times branching ratio a{pp —)• tH^ -|- X) x BR(Ff^ —has been given to be between 


0.76 pb and 4.5 fb in the range of mH± =180 GeV to 1 TeV 


36l |. This limit gives an upper limit 


on tan/3 for a fixed value of mij± in the 2HDMs. For example, tan/3 > 50 (60) at mu± = 200 
(230) GeV can be excluded at 95% CL in the MSSM [36|], where a similar bound is expected to 
be obtained in the Type-II 2HDM because of the same structure of the Yukawa interaction®. In 
the Type-I and Type-X 2HDMs, the production cross section of pp —)• -|- V is significantly 

suppressed by a factor cot^ /3, so that we cannot expect to obtain an important bound in the high 
mass region. 


G. Summary of the constraints on 

In Tab. m we present the summary of the current experimental bounds on mu± in the 
2HDMs/3HDMs with the four types of Yukawa interactions from various experimental observa¬ 
tions. 


V. NUMERICAL RESULTS 


In this section, we perform numerical evaluations for the H^W^V vertices and related ob¬ 
servables. In particular, we focus on the light case, i.e, mfj± = 0(100) GeV, because of its 
phenomenological interest. As we discussed in Sec. IV, such a scenario is allowed in the Type-I 
and Type-X Yukawa interactions from flavour constraints, so that we consider these types only in 


^ Notice that we have emulated both the t op quark production and the decay as well as R^-strahlung through the 
single gg —>■ tbH^ mode, in the spirit of 1371] . 

® In the Type-Y 2HDM, although the same production cross section of pp —>■ -|- Y is obtained as in the Type-II 

case, the branching fraction of —>■ is significantly suppressed due to the enhancement of the decay rate of 
the —>■ be mode [s^. Therefore, the bound in the Type-Y 2HDM can be much weaker than that in the Type-II 

case. 
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Experiment 

95% CL lower lim. on mH± 

tan/3 

Type 

Comments 

6 —?► 57 

322 GeV 

- 

II and Y 



(800, 200, 100) GeV 

(1, 2, 2.5) 

I and X 


B°-B° 

(500, 300, 100) GeV 

(1, 1.5, 2) 

All 


LEP II 

(80, 90) GeV 

- 

All 

Brv + Beg = 1, Bti/ = 1 

t H^b 

(160, 140, 100) GeV 

(L 2 , 4) 

I 

Using 1.3% (See Fig.[2|) 

at the LHC Run-I 

(160, 150, 130) GeV 

(L 2 , 4) 

X 

Using 1.3% (See Fig.[2|) 


TABLE II: The 95% CL lower bound on mH± in the 2HDMs/3HDMs from various experimental mea¬ 
surements for a fixed value of tan/3. For the row of LEP II, 80 (90) GeV is given for the case of 
Bru + Bcs = 1, {Btv = 1), where Bru and Bcs are the branching fractions of and ^ cs 

modes, respectively. 


this section. First, we evaluate the form factors of the H^W^Z and vertices. For the 

vertex, since the form factor is related to G-^ by the Ward identity, we only show 
and Second, we show all the branching fractions of , including the —?> W^Z and 

7 modes. Finally, we discuss cross sections for various signal processes involving the 
^W^Z and decays at the LHC. 

In our model, there are 16 independent parameters in the potential given in Eq. ([T]), namely, 
/i^, Ai_ 5 , Xri, / 3 i -3 and (Ti_ 3 . They are divided into 8 parameters in the active sector (/x ^_3 and 
Ai_ 5 ) and the remaining 8 parameters (/i^, Xri, pi-z and ( 7 i_ 3 ). 

After the tadpole conditions are imposed, the former 8 parameters can be expressed by v, tan /3, 
sin(/3 — a) nih, rrijj, m^, mjj± and M^. Two of the 8 parameters, v and nih, should be used to 
reproduce the gauge boson masses and the observed Higgs boson mass, i.e., v ~ 246 GeV and 
nih — 125 GeV. Furthermore, the Higgs boson search data at the LHC suggests that the observed 
Higgs boson is SM-like [l-4|, so that taking sin(/3 — a) ~ I gives a good benchmark scenario as we 
explained in Sec. H. We thus take sin(/3 — a) = 1 in the following calculation. 

Regarding the latter 8 parameters, we proceed as follows. First, we take Xrj = 0, as this gives 
a four-point interaction among the inert scalar bosons that does not affect the following analysis. 
Second, we take pi and ui so as to satisfy the vacuum stability condition given in Eqs. (1411) and 
([4^ for given values of p 2,3 and < 72 , 3 : 


Pi = MIN( 0,/92 + P3,P2 - P 3 ), cTi = MIN( 0 ,U 2 -b (73, a 2 - 0 - 3 ). 


(61) 
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Finally, the remaining 5 parameters can be expressed in terms of three masses of the inert scalar 
bosons {m^±, rrirf^ and and the p 2 and ps parameters. In this parametrisation, the <72 and 

(73 parameters are given as the outputs: 


<72 = -P2 cot^ P + 


0-3 = -P3 cor P + 


1 


sin^ P 

1 

sin^ P 


- 2m2±) , 


(62) 

(63) 


Therefore, to recap, we are left with 5 new parameters in the active sector (tan/3, m//±, 
rrijj and M^) and 5 new ones in the inert sector too (m^±, rn„ , ^ and ps) and we will scan 


over these. Regarding the SM inputs, we use the following values 


39|, 


40]: 


mz = 91.1876 GeV, mw = 80.385 GeV, Gp = 1.1663787 x 10"^ GeV^ 

mt = 173.07 GeV, rub = 3.0 GeV, ruc = 0.677 GeV, Vcb = 0.0409, Vts = 0.0429, 

rur = 1.77684 GeV, = 0.105658367 GeV, nih = 125 GeV. (64) 


where Vij are the Gabibbo-Kobayashi-Maskawa matrix elements, and the quark masses mb and me 
are given at the m^ scale as quoted from Ref. Q- 

The form factors depend on the three momenta p^, py and + Py ^ (= -^jT) 

and respectively. In the numerical calculation, when mfj± > myr + m 2 ;, we take p^ = m^, 
p ^2 = ~ while when mu± < m^r + we take p^ = {mp± — mz)^, p% ~ 

and (thereby allowing for below threshold decays too). For the vertex, we 

take p^ = = 0 and q^ = ut,^± . 


A. Form factors 

We start by showing the numerical results of the form factors of the H^W^Z and 
vertices. In order to see how the inert scalar boson loops can change the prediction, we first 
show the result in the 2HDM under the constraints from unitarity, vacuum stability and the EW 
parameters as discussed in Sec. IV. Then, we move on to the 3HDM. 

In Fig. [3l the values of \Xz\^ {X = F, G and H) and \Yp^ {Y = G and H) are respectively 
plotted in the upper and lower panels as a function of tan/3 in the case of sin(/3 — a) = 1 and 
mjj = mH±. The left (right) panel shows the case of mH± = 150 (200) GeV. The solid, dashed 
and dotted (dashed and dotted) curves respectively show the fermion loop contribution to \Fz\^, 
\Gz\^ and \Hz\‘^ (I^^P and \HP^) whereas the black and blue (blue) scatter plots are the total 
contribution to \Fz\^ and \Gz\^ (|G..j.p), respectively. For the boson loop contribution, we scan the 
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FIG. 3: Values of {X = F, G and F[) (upper panels) and \Y^\^ {Y = G and H) (lower panels) as 

a function of tan/3 in the 2HDM with the Type-I or Type-X Yukawa interactions. We take = 150 
GeV (left panels) and 200 GeV (right panels). In both the panels, mn = 'rnu± and sin(/3 — a) = 1 
are taken. The values of and m\ are scanned over the ranges of —400^ < < +400^ GeV^ and 

100 < niA < 260 (350) GeV in the left (right) panels, respectively. The solid, dashed and dotted (dashed 
and dotted) curves respectively show the fermion loop contribution to \Fz\^, |Gzp and (IG^P and 

while the scatter plots show the total contribution. 


parameters over the intervals —400^ GeV^ < < 400^ GeV^ and 100 GeV < < 260 (350) 

GeV in the left (right) panels. We note that > 260 (350) GeV when mjj±{= nij^) = 150 (200) 
GeV is excluded by the constraint from the S parameter at 95% GL. We also note that only the 
fermion loop contributes to Hz and 

We can see that the value of \Fz\'^ is the biggest of all the form factors as we expected in 
Sec. Ill, because of the dependence. Typically, \Fz\'^ is more than one order of magnitude 
larger than \Gz\'^ and \Hz\‘^- In addition, all the squared form factors decrease as tan/3 is getting 
larger, because the top Yukawa coupling is proportional to cot /3. The maximal allowed value of 
is obtained to be about 10“^ at tan/3 ~ 2.5 in both the cases of mfj± = 150 GeV and 200 
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FIG. 4: Values of (X = F, G and H) (upper panels) and |V-yp {Y = G and H) (lower panels) as 

a function of tan /? in the 3HDM with = 400 GeV. We take = 150 GeV (left panels) and 200 
GeV (right panels). All the other parameters are taken as given in Eq. (1551) . In the upper panel, the black 
scatter plot shows the values of |Fzp. In all the panels, the blue scatter plot and the solid curve respectively 
represent |Gyp and {V = Z,^). 


GeV. For the vertex, the maximal allowed values of and \H^\^ are order of 10“® at 

tan /3 ~ 2. 

Regarding the 3HDM, as we see from Eq. (|37p . Fz is logarithmically enhanced by in the 
case of m^± = However, a too large mass difference between and t]^ is excluded by the 

S parameter as shown in Fig. [T]in the case of mjj± = = rUfj or AT = 0. We thus take a mass 

difference between F[^ and A/H with to avoid the constraint by the effect of non-zero 

AT. From the above reason, we consider the following parameter conditions in the forthcoming 
calculations: 


= rrijj = mjj± + 50 GeV, , 

1 , , 
= TUr,^ = -m^, TUr,^ > m^±, -10 < P2,P3 < 10. (65) 
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FIG. 5: Values of {X = F, G and F[) (upper panels) and \Y^\^ (Y = G and iJ) (lower panels) as a 

function of in the 3HDM with tan/3 = 2.5. We take = 150 GeV (left panels) and 200 GeV (right 
panels). All the other parameters are taken as given in Eq. (1551) . In the upper panel, the black scatter plot 
shows the values of |Fzp. In all the panels, the blue scatter plot and the solid curve respectively represent 
|Gy|2 and \Hv\^ {V = Z,^). 


We note that, in this setnp, r}^ corresponds to the DM candidate. The measured relic abundance 
of DM® can be satisfied by the resonant process of rj^r]^ —)■ A/H —//. 

In Fig. m the values of \Xz\‘^ {X = F, G and H) and |KyP {Y = G and H) are respectively 
shown in the upper and lower panels as a function of tan /3 with = 400 GeV. The left (right) 
panel shows the case of mjj± = 150 (200) GeV. In the upper panel, the black scatter plots show the 
values of iF^p. In all the panels, the blue scatter plot and the solid curve respectively represent 
\Gv\‘^ and \Hy\'^ {V = Z,^). Similar to the results in the 2HDM, \Fz\^ is the biggest of all the 
squared form factors also in the 3HDM, and all the squared form factors become smaller when 
tan/3 becomes large. Remarkably, at tan/3 = 2, we obtain \Fz\^ — 10“^, which is one order of 

® Because the DM phenomenology is not the main topic of this paper, we do not perform the detailed analysis such 
as the calculation of the (co)annihilation cross sections of the DM candidate. 
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magnitude larger than \Fz'^ in the 2HDM. 

In Fig. O we show the dependence of the squared form factors in the case of tan [3 = 2.5. 
We take = 150 (200) GeV in the left (right) panel. The description of the objects in the 

figure is the same as in Fig. 01 Clearly, we can see that only \Fz\‘^ is enhanced as is getting 
larger. The maximal allowed value of \Fz\‘^ is about 10“^ at ~ 500 GeV. 


B. Branching fractions of 


Next, we discuss the decay branching ratios of H^. As we see in Figs. 0] and [5] that the form 
factor Fz is much larger than Gz and Hz, we only keep the term proportional to \Fz\‘^ for the 
WZ decay. When > mw + rnz, the on-shell decay of W^Z opens and its 

decay rate is calculated as 

—)■ W^Z) = + ^‘^xy^xz]\Fz\^, ( 66 ) 

where ~ If smaller than the off-shell decay 

modes —)■ W^Z* and —>• W^*Z are allowed. The decay rate with three body final states 

is given by 


^W^Z* ^ Wff) 
f 


2567r3mj|^± 


\Fz?‘F^ {Xz,Xyr) , 


^g‘^rn\ 

bl2'K'^m}j± 


\Fz? 



40 9 



F^ {Xyr,X2:) , 


(67) 

( 68 ) 


where 


Fz{x,y*) 

1 

2Axy* 


arctan 


y* {l+x) — {l—xY 


4xy^—A(x, y*) 
{(x-l)[6y*2+y*(39x- 


T 71" 

-[(1 - y*f - 3x3 ^ ^ 

9) + 2(l-x)2]-3y*[y*' + 22/*(3x- 


-^{l-y*fx 

1) — x(3x -b 4) -|- 1] Inxj. 


(69) 


We note that the argument y* is for the ratio of squared masses of a virtual gauge boson to that 
of e.g., for the —)■ W^*Z case, we should use F^{m\/m^^±,rn^/rn?^±). The decay rate 

for 7 is given by 

r(i7± ^ w±7) = - ^ty)' {\G,\^ + ■ (70) 

In Fig. [6l we show the branching fractions of as a function of in the 3HDM with the 
Type-I Yukawa interaction. We take mu± = 150 (left), 170 (center) and 200 GeV (right). The 
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FIG. 6: Branching fractions of as a function of in the Type-I Yukawa interaction with tan j3 = 2.5. 
We take = 150 GeV (left), 170 GeV (center) and 200 GeV (right). 


+ 

S’ 

CQ 



™I1A ™I1A ™t1A 

FIG. 7: Branching fractions of as a function of in the Type-X Yukawa interaction with tan/3 = 2.5. 
We take = 150 GeV (left), 170 GeV (center) and 200 GeV (right). 


value of tan /3 is fixed to be 2.5 in all the panels. In these plots, we scan the values of p 2 and ps in 
the range of —10 to +10 and extract the set of {p 2 ,P 3 ) combinations giving the maximal value of 
the decay rate r(i7^ —)• WZ). Further, for the case of < mw + rnz^ we show the branching 
fraction of —)• W^Z as the sum of the branching fractions of —)• W^Z* and —)• W^*Z. 

In all the plots, the behavior of in the > W^Z decay is similar to that of \Fz\^ shown in 
Fig. O In the case of = 150 GeV, although BR(FI^ —^ W^Z) benefits from the enhancement 
of |i^zp, its rate is smaller than BR(FI^ —)• VF^y) when < 300 GeV. This can be understood 
by the suppression of the decay rate of ^ W^Z due to the off-shell effect of the or Z 

bosons. Therefore, we obtain a larger value of BR(FI^ —>• W^Z) in the case of mfj± = 170 GeV 
because of the smaller off-shell effect. However, once mij± exceeds the top quark mass, both the 
branching fractions of —>• W^Z and —>• FF^y are significantly suppressed by the —)■ tb 

decay. We find that the maximal value of BR{H^ —)• W^Z) is about 4%, 40% and 0.4% in the 
cases of m}j± = 150, 170 and 200 GeV, respectively. 

In Fig. [71 we also show the branching fraction of in the Type-X Yukawa interaction with 
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Type-I 

Type-X 

Br(t ^ HH) [%] 

( 3.3, I.IO, 4.7x10- 

-3) (3.3, l.l, 4.7x10-3) 

Br(iJ± ^ W^Z) [%] 

(0.66, 3.5, 33) 

(0.025, 0.14, 1.8 ) 

Br(i7± ^ 1T±7) [%] 

(1.6, 2.1, 1.6) 

(0.059, 0.081, 0.087) 


(390, 700, 29) 

(15, 28, 1.6) 

<^5:^ [fb] 

(940, 420, 1.4) 

(35, 16, 0.075) 

[fb] 

(2.3, 7.5, 46) 

(0.087, 0.30, 2.5) 

[fb] 

(5.5, 4.5, 2.2) 

(0.20, 0.17, 0.12) 


TABLE III: The branching fractions and the cross sections in the 3HDM with a Type-I and Type-X Yukawa 
interaction. We take tan/3 = 2.5 and = 400 GeV. The numbers in the bracket correspond to the result 
of =130, 150 and 170 GeV from left to right. 



Type-I 

Type-X 

Br(t ^ H^b) [%] 

(1.3, 0.43, 1.8x10- 

-3) (1.3, 0.43, 1.8x10-3) 

Br(i7± ^ W^Z) [%] 

(0.52, 2.7, 26) 

(3.0x10-3, 0.016, 0.21) 

Br(H± ^ 1T±7) [%] 

(l.l, 1.5, 1.2) 

(6.5x10-3, 8.6x10-3, 9.3 x10-3) 

[fb] 

(120, 210, 8.6) 

(0.71, 1.3, 0.070) 

[fb] 

(260, 120, 0.40) 

(1.5, 0.68, 3.1x10-3) 

J [fb] 

(1.8, 5.8, 36) 

(O.OlO, 0.034, 0.29) 

[fb] 

(3.8, 3.2, 1.7) 

(0.022, 0.018, 0.013) 


TABLE IV: Same as Table ITTTl but for tan/3 = 4. 


tan/3 = 2.5. Although we observe a similar behavior of —)• W^Z) and —)■ 

as seen in Fig. [6l their maximal values are smaller than those in the case of the Type-I Yukawa 
interaction. This is because in the Type-X Yukawa interaction, the decay rate of the —>• 
mode is enhanced by tan^ /3. Here, the maximal value of —>• W^Z) is about 0.2%, 2% and 

0.3% in the cases of = 150, 170 and 200 GeV, respectively. 
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C. Cross sections at the LHC 


Finally, we discuss the signature of the —)• W^Z and —)■ decays at the LHC. If the 

mass is below the top quark mass, the top decay t —)• H^b is the dominant production mode 
of while above it i/^-strahlung becomes dominant. In reality, the latter is never significant 
as a means of enabling —)• W^Z and —)• IF^y detection, so we only concentrate on the 

former. We then expect the signature pp —?■ bbH^W^ —)• bbW^W^V. The signal cross section of 
this process is estimated by 


fjgP = 2 X ati X [1 - BR(t ^ HH)] x BR(t ^ H^b) x BR{H^ W^V), 


r± 


±1 


(71) 


where a^i is the top quark pair production cross section at the LHC. In Ref. 


■ Q. 


ati = 923.0 pb 


has been obtained with rrit = 171 GeV and = 14 TeV at the next-to-next-to leading order 


using CTEQ6.6 parton distribution function 


421. As alternative production modes of states. 


especially helpful when the charged Higgs mass is larger than the top quark mass, one should also 
count the EW productions, e.g., pp H^A, pp —)• and pp —)• whose cross sections are 

determined by the masses of extra Higgs bosons. The cross sections for H^A and productions 

are the same as long as we take and sin(/3 — a) = 1. By using these production modes, 

we can consider pp — ^ H^A/H^H — W^V + and pp — >• — )■ W^V + A^, where 

and X^ are respectively the decay product of AjH and 77^. The signal cross section via the EW 
production modes are estimated by 


+ 2aH+H-) X BR(i7± ^ IF±R), (72) 

where (Th^a^ cth+h- S'I'® respectively the cross sections of pp -A- H^A, pp -A and 

pp -A 77+77“. In the cases of m}j± = 130, 150 and 170 GeV, we obtain (Jh±a {(^h+h-) = 84 (89), 
54 (53) and 36 (34) fb , respectively, at y/s = 14 TeV using CTEQ6L. Eor crjj±j^ (= crjj±jj), the 
above numbers are obtained by summing the H^A and H~ A processes. 

In Tabs, [m] and lYl we show the branching fractions of the t -A H^b, 77+ -A W^Z and 
77+ -A IF +7 modes and the overall signal cross sections of both the top decay and EW processes 
estimated by using Eqs. dm) and (El]), respectively. The results with the Type-I (X) Yukawa 
interaction are given in Tab. IHII (lIVp . For the top decay process, the production cross section 
gets smaller when m[j± approaches nit because of the phase space suppression. Conversely, the 
branching fraction for 77+ -A VF+Z becomes larger as we already seen in Figs. [6] and [71 As a 
result, (7^°^ attains a maximal value around ^150 GeV, while is simply reduced as 
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m^± becomes larger since BK{H^ —)• does not encounter any threshold (as mjj± > m\Y±). 

For the EW processes, the reduction of the production cross section {(7h±aj ^h±h ^h+h-) is 
milder than that of the top decay process ((TjfxBr(t —?■ H^b)). Therefore, the signal cross section 
of the EW processes become larger than the top decay process at mH± = 170 GeV. Einally, we 
note that the signal cross sections in the Type-X case is more than one order of magnitude smaller 
than those in the Type-I case. 

VI. CONCLUSION 

We have computed the strength of the H^W^Z and vertices at the one-loop level in 

the 3HDM under a Z 2 x Z 2 symmetry, which defines a Higgs sector with two active doublets and 
one inert one. We have discussed all the four types of the Yukawa interactions which are defined by 
the Z 2 charge assignment to the SM fermions. We have taken into account vacuum stability and 
perturbative unitarity as theoretical constraints, and have considered the bounds from the EW 5, 
T and U parameters, flavour experiments and direct searches for states at LEP-II and LHC 
Run-I. We have seen that the mass of the can be smaller than the top quark mass in models 
with the Type-I and Type-X Yukawa interactions, but not in Type-H and Type-Y. Eurther, we 
have shown that, among all the form factors, only Fz can be enhanced with respect to the 2HDM 
by taking large mass splittings between and 7 /^, because of the non-decoupling effect of the 
inert scalar boson loop contributions. 

In particular, we have found that in the 3HDM the squared form factor |Ezp can be one 
order of magnitude larger than that predicted in the 2HDM under the aforementioned theoretical 
and experimental constraints. In addition, the branching fraction of the ^ W^Z mode 

can be about 4 ( 0 . 2 )%, 40 ( 2 )% and 0.4 (0.3)% in the cases of mu± = 150, 170 and 200 GeV, 
respectively with the Type-I (Type-X) Yukawa interactions. In contrast, the branching fraction of 
the 7 mode is at the few percent level as long as mff± is smaller than the top quark 

mass in the Type-I and Type-X cases, thus benefiting from very little enhancement with respect 
to the 2HDM. Such increased rates in the 3HDM stem from loop contributions due to inert Higgs 
states that are absent in the 2HDM. 

Finally, we have discussed signal processes embedding —)• W^Z and —)• decays 

at the LHG. In the light scenario, i.e., mu± < mt, with the Type-I and Type-X Yukawa 
interactions, the top quark decay process t —)• H^b is the dominant production mode for 
except for the extreme case of Frut- In the heavy scenario, i.e., mj{± > mt, this channel 


26 


is no longer viable and we have resorted to the bg —)• tH^ mode. (Herein, we have emulated top 
production plus decay and ff=*=-strahlung via gg —)• tbH^.) In fact, there are also EW production 
modes, such as pp —pp —)• H^A and pp —?■ By combining the production and 

decay of we have considered the signal processes pp —)• bbW^H^ —)• bbW~^W~V, pp —>■ 

H~^H~ —^ W^VX^ and pp —)• H^A/H^H —)• W^VX^. We have thus computed the ensuing 
cross sections in all cases and shown that the LHC Run-H has the potential to access —)• W^Z 

and/or —)• decays, certainly for light LI^’s (at standard luminosity) and possibly for 

heavy LI^’s (at very high luminosity). To establish one or the other such signals at the CERN 
machine may represent circumstantial evidence of a 3HDM sector, as opposed to a 2HDM. 
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Appendix A: IPI contributions 


Here, we give the analytic expressions for the IPI diagram contribntions to the form factors 
of the H^W^V iy = Z, 7 ) vertices and those for the W^-H^ mixing The fermion 

loop contribntion to the H^W^V vertices has been calcnlated in Ref. 1C I whereas the boson 

□ □ 

contribution in the 2HDM has been evaluated in Refs. [111.Il2l|. In addition to these contributions, 
there are inert scalar boson loop contributions as shown in Fig. [HI 

In the following, we separately show the fermion and boson loop contributions to the form 
factors denoted by X^p and X^p (A = F,G and H), respectively. Regarding the boson loop 
contribution, we only show the contributions from pure scalar loop diagrams, where scalar bosons 
are running in the triangle and circle type diagrams (see Fig. [ 8 ]). There are additional gauge-scalar 
mixed type diagrams, where one gauge and two active scalar bosons or two gauge and one active 
scalar bosons run in the triangle part. Because these contributions are proportional to cos(/3 — a), 
they vanish or become negligible by taking the SM-like limit sin(/? — a) —)• 1 or taking the SM-like 
regime sin(/3 — a) ~ 1, respectively. We thus neglect them here^^. 
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FIG. 8: The IPI diagrams for the H^W^Z and vertices. The diagrams which vanish in the limit 

sin(/3 — a) = 1 are not displayed. 


The contributions from the gauge-scalar mixed type diagrams are given in Ref. [n 
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FIG. 9: Diagrams giving the H^-W^ mixing. The diagrams which vanish in the limit sin(/3 — a) = 1 are 
not displayed. 






In order to express loop functions, we use the Passarino-Veltman functions [43|]. Here, we give 
the integral formulae of some of the functions which we use in the following discussion: 

Ho(p^;mi,m2) = A - / (ixlnA^, 

Jo 

A 

Bi{p‘^;mi,m 2 ) = -— + / iix(l - x) In 
^ Jo 

Co{pl,phq^]'mi,m 2 ,m‘i) = - I dx dy^, 

Jo Jo Ac 

2 2 2 idif-ii 

Cu{Pi,P2^Q ;mi,m2,m3) = - / dx dy 


Ac ’ 
y{y - 1 ) 

Ac ’ 

2/(1 - xy)"^ 

Ac ’ 

. 2/( 1 - yf 

Ac 

2/(1 - xy){l - y) 


c 


A 1 /■! /•! 

C'24(p?,pi,g^;mi,m2,m3) = — --y dx j dyyln 


Ac, 


(Ala) 

(Alb) 

(Ale) 

(Aid) 

(Ale) 

(Alf) 

(Alg) 

(Alh) 

(Ali) 


where 


As = —x(l — x)p^ + xm\ + (1 — x)m 2 , 

Ac = y'^iPix +P2f +y[x{p2 -q^ + ml- ml) + ml - ml - pi] + ml, 


(A2) 

(A3) 


In Eq. ()Alh . A is given by 


A=- 7 c + In dvr + In(A4) 

e 

where e appears in the D{= 4 — 2e) dimensional integral, p is an arbitrary dimensionful parameter 
and 7 ^ is the Euler constant. In the four dimension limit e —^ 0, A is divergent. We note 
that this divergent part A appears in the following expressions, but it is exactly cancelled in the 
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renormalized variables such as Xz and {X = F, G and H). We use the shorthand notations 


The fermion loop contribution to X^^ is given by 


iPi _ 2Nc 


like Bi{p^-,A,B) = Bi{p^-,mA,mB) and Cj, ij{A,B,C) = Ci^ ij{pf,p^,q^-,mA,mB,mc). 

IP 

z 

w ^ 

+ rnf^tivb + Oft) 4C24(t, b, b) - Bo{q‘^; mt, mb) - Bo{pl^-mb, mt) - {2ml “ Pz)Co{t, b, b) 

- mlib{vb + Ob) 4C'24(t, b, b) - Bo{p%] mb, mb) - Bo{q‘^-, mt, mb) - {ml + ml - p^)C'o(t, b, b) 

- mlib{vb - ab) BQ{p%;mb, mb) + Bo{ply;mt,mb) + {ml + ml - q'^)Co{t, b, b) 

+ 2mlmlit{vb - Ob)C'o(t,6,5)| + {mt,it,Vb,ab) ^ {mb,-(b,vt,at), 

= 2 2^ '>^t^t{vb + Oft)(2(723 + 2(7 i2 + Cii + C(j) 

iovr V L 


(A5) 


-m^^ft(uft +Oft)(2(723 +C 12 ) - mlCb{vb - ab){Ci2 - Cn) {t,b,b) 
+ {mt,Ct,Vb, ab) ^ {mb, -^b, Vt, at). 


(A6) 


^ 4iVeWt^ 


rrlPI _ 
ri 7 T? — 


W 


mlit{vb + aft)((7o + Cu) - mlib{vb + Oft)(7i2 + mlib{vb - aft)((7i2 - Cn) {t, b, b) 


+ {mt,^t, Vb, ab) o {mb, +6, vt, at). 


where 


W = 2^f - ^wQf^ af = -If. 


(A7) 


(A8) 


That to X}f^ is given by 


‘^XcQb 


= 


w 


{ 


+ ml^t 4C24(t, b, b) - Bo{q‘^;mt, mb) - -Bo(p^; mb, mt) - (2mf - p^)(7o(t, b, b) 

- mlib 4(724(t, b, b) - Bq{p\-, mb, mb) - Bo{q‘^;mt, mb) - {ml + ml - Pw)Co{t, b, b) 

- mlib Bo{p\-,mb, mb) + Fo(p^; mt, mb) + {ml + ml - q^)Co{t, b, b) 

+ ‘^rnlmlCtCo{t, b, 6)| + {mt, Qb) ^ {mb, -^b, Qt), 


(A9) 


Gi^l = 


2 r 


ANcQbm^^r 
IGvr^u^c 


w 


'Oi??t(2(723 + 2 Ci2 + Gii + Cq) 


- "i6Cb(2(723 + G 12 ) - mlib{Ci2 - Gii) {t, b, b) + {mt,Ct, Qb) ^ {mb, -?ft, Qt), (AlO) 
2 r 


TjiPi _ ^^cQbrniw 


mlit{Gt) + (7ii) - mlibCi2 + mlib{Ci2 - Gn) {t, b, b) 
+ {mt, 6, Qb) ^ {mb, +^b, Qt)- 


(All) 
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The boson loop contribution is given by 


77 IPI _ ^ f 

^z.B — TTT-^ - ) 

T V 


VoTx'^VC^ 

+ ^mu-H sm(/3 - «) [(2 - 4 s^)C24(^, - 2 C 24 {H^,A, H) + sl,Bo{q^-, H^,H) 
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(2 - As'^)C24ivA,V^,V^) - ‘^C24iv^,riH,ri^) + s‘^Bo{q^]ri^,rij^) |, 
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“ ^H+ri- 


Va 


(A12) 


/-»iPi _ 
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IbTT^UC^ 

+ ^H+H-H sm(/3 - a) [(2 - 4 s^)(C 7 i 2 + C 23 ){H, H^,H^)- 2{Ci2 + C 23 ){H^,A, H) 

(2 - 4s^)(C'i2 + C23){'nH,V'^,V^) - 2(^12 + C23){v^,Va^ Vh) 

(2 - 4s^)(C'i2 + C23){r]A^^f ^r]^) “ 2(^12 + C23){v^,VH,r]A) }, 


- ^H+ri-r,„ 


-A 


^H+rj-r/^ 


(A13) 


pipi _ 
^■y.B — 


ATs = - a)[AC24iH,H^,H^) - Bo{q‘^-H^,H)] 

[46*24(??H, Bo{q‘^;iq^, rj^)] - A^+r,-,,^ [4C24{Va^ Bo{q‘^] r/=^, r/^)]}, 

(A14) 
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“ ^H+rj- 


TTT^ sm (/3 - a)(Ci2 + C23){H, 

lOTT^V L 

(6*12 + 6*23)(r?//,r/=^,r/=^) - A^^+^-^^(Ci2 + C'23)(r?^, r/=^, ??=^) 


■nH 
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and 


rrlPI 

^Z,B 


= = 0 


(A16) 


where 


^H+H-H 


- (rnjr — M^)(cot /3 — tan / 3 ) sin (/3 — a) — { 2 m\± + 

V \- 

V 

-^{P2 + /03 - 0-2 - 0-3) sin 2 / 3 , 

V 

±^(P2 -P3 -0-2+ 0-3) sin 2 / 3 . 


2 M^) cos(/3 



(A17) 


(A18) 
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We note that the above expressions are obtained by extracting the coefficient of the scalar trilinear 
vertex, i.e., £ = 4 -• 

The fermion and boson loop contributions to the W^-H^ mixing, i.e., and 
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^wh{P)Bi respectively, are given by: 

^wh{p^)f = + Bi) - mlibB]\{p‘^]mt,mb), 


Ibvr^ 

rlPI (r?\ — * '^W 

r„-„(p )b - r^ — 


A/i-+/f-Hsin(/3 - a){2Bi + BQ){p^;mH±,rnH) 


+ + *Vr;-r?A( 2 ^i + ^o)(p^;%±,m^^) 

The counter term contribution is then obtained from the above W^-H^ mixing via Eq. 
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4s^-/Vc r™2t /D I D ^ D l/„2. 


lOvr^u^c 




w 


’W 


Ibvr^ 


uc 


W 
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SR 
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(A20) 

(A 21 ) 

(A 22 ) 

(A23) 

(A24) 


Using the above analytic expressions, we can directly check the relation from the Ward identity 
in Eq. ([32]), i.e., + 5F^) = - m^/m^±)/2. 
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